Sequential tunneling and inelastic cotunneling in nanoparticle arrays 
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We investigate transport in weakly-coupled metal nanoparticle arrays, focusing on the regime 
where tunneling is competing with strong single electron charging effects. This competition gives 
rise to an interplay between two types of charge transport. In sequential tunneling, transport is 
dominated by independent electron hops from a particle to its nearest neighbor along the current 
path. In inelastic cotunneling, transport is dominated by cooperative, multi-electron hops that 
each go to the nearest neighbor but are synchronized to move charge over distances of several 
particles. In order to test how the temperature-dependent cotunnel distance affects the current- 
voltage (J — V) characteristics we perform a series of systematic experiments on highly-ordered, 
close-packed nanoparticle arrays. The arrays consist of ~ 5.5nm diameter gold nanocrystals with 
tight size dispersion, spaced ~ 1.7nm apart by interdigitating shells of dodecanethiol ligands. We 
present I — V data for mono-, bi-, tri- and tetralayer arrays. For stacks 2-4 layers thick we compare 
in-plane measurements with data for vertical transport, perpendicular to the array plane. Our 
results support a picture whereby transport inside the Coulomb blockade regime occurs by inelastic 
cotunneling, while sequential tunneling takes over at large bias above the global Coulomb blockade 
threshold Vt(T), and at high temperatures. For the smallest measurable voltages, our data is fit 
well by recent predictions for the temperature dependence zero-bias conductance due to multiple 
cotunneling. At finite, but small bias, the cotunnel distance is predicted to set the curvature of the 
nonlinear I—V characteristics, in good agreement with our data. The absence of significant magnetic 
field dependence up to 10T in the measured I — V characteristics further supports the picture of 
inelastic cotunneling events where individual electrons hop no further than the nearest neighbor. 
At large bias, above the global Coulomb blockade threshold, the I — V characteristics follow power 
law behavior with temperature-independent exponent close to 2 predicted for sequential tunneling 
along branching paths that optimize the overall charging energy cost. 



I. INTRODUCTION 

Arrays of closely-spaced nanoparticles are currently be- 
coming a model system for fundamental studies of meso- 
scopic charge transport and for targeted applications in 
nanotcchnologyi Quite generally, their behavior depends 
on the competition of several energy scales that can be 
controlled independently by tuning the properties of the 
particles and their geometric arrangement. These are the 
mean energy level spacing, the Coulomb charging energy 
for a single particle, as well as the tunnel energy associ- 
ated with the inter-particle coupling^ 

Even for particles that individually are simple met- 
als, the array as a whole can exhibit the full range 
from strong, exponential ("insulating") to weak, non- 
exponential (" metallic" ) temperature dependence of the 
conductance. This behavior can be tuned by the strength 
of coupling between the building blocks - the bare, high- 
temperature tunneling conductance g. In the regime 
of strong coupling, the Coulomb interaction is screened 
and electrons propagate easily. In the opposite limit 
of the weak coupling, single electron charging becomes 
significant, leading to Coulomb blockade behavior and 
the localization of electrons on individual particles^ As 
a results, arrays with strong coupling g > g q , where 
g q = e 2 /h is the quantum conductance, behave as met- 
als, while arrays with g < g q are insulators in the low- 



temperature, low-bias limit. 

The overall, global transport properties aris- 
ing from local, charging-mediated tunneling have 
been investigated experimentally for a number of 
systems formed from metallic or semiconducting 
nanoparticles.i^^^^^^&i 7 -^ Close-packed 
two-dimensional arrays of highly size-controlled metal 
nanocrystals, in particular, have served to elucidate the 
effect of structural disorder, such as imperfections in the 
particle arrangements, on the current-voltage character- 
istics^i 1 ^ and to test predictions 33 for quenched charge 
disorder i^ii Charge disorder arises from variations in the 
local chemical potentials, due to polarization by trapped 
parasitic charges in the substrate or in the ligand shells 
surrounding the particles. In practice, charge disorder 
is unavoidable for arrays involving more than a few 
particles and, as a consequence, tunneling occurs along 
paths that optimize the overall energy cost. This results 
in two distinct regimes. At applied bias voltages large 
enough to overcome local Coulomb energy costs, trans- 
port occurs by sequential tunneling between neighboring 
particles along a set of optimal paths. 17 ! 33 At small bias 
and low temperatures, sequential tunneling is suppressed 
by the Coulomb blockade. In this regime, conduction 
involves higher-order, cooperative tunneling processes, 
so-called cotunneling events, that can transport charge 
over distances of several particles without incurring 
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the full Coulomb energy costsi 2 ' 27 ' 28 ' 29 By virtue of 
the uncertainty principle, this is possible through the 
creation of short-lived excitations on the intermediate 
particles. In arrays of metallic nanoparticles and for 
temperatures from a few Kelvin on up, these excitations 
consist of electron-hole pairs, produced by cotunncling 
events in which one electron tunnels onto a particular 
particle while simultaneously another, lower energy 
electron tunnels off (inelastic cotunneling) . 

Recently we investigated electron transport through 
large, highly-ordered metal nanoparticle arrays with lat- 
eral size of 65-70 particles between the electrodes.- 2 We 
found that the resulting, current-voltage characteristics 
as well as the temperature dependence of the conduc- 
tance for small applied bias voltages are in line with pre- 
dictions based on multiple inelastic cotunnelingi^ i 27 ! 29 In 
particular, our results indicated that low-temperature, 
low-bias transport across the full array occurs via a se- 
quence of cotunneling events, each involving a few par- 
ticles (up to ~ 4), depending on temperature. This 
temperature-dependent cotunneling distance was inter- 
preted as resulting from an optimization process that 
weighs the likelihood of n-electron, cooperative cotun- 
neling events against the net energy cost. 

In the present paper we extend these results and specif- 
ically address the question about what sets the cotunncl- 
ing distance. This characteristic length scale not only de- 
termines the low-bias current-voltage characteristics, but 
also the temperature dependence of the zero-bias conduc- 
tance. We present results from experiments in which the 
distance between electrodes has been tuned to become 
commensurate with the cotunneling distance. These ex- 
periments involve vertical transport through carefully 
prepared stacks of particle layers, up to four layers in 
height. This approach provides a means to control the 
number of tunnel junctions the electron traverses to get 
across the electrodes, allowing us to test the cotunneling 
scenario in detail. 

This paper is organized as follows. In Sec.|TT]we discuss 
general features of the current- voltage (I — V) character- 
istics predicted for metal nanoparticle arrays. Sec. IIIII 
contains details about the sample fabrication and the ex- 
perimental set-up. We present in Sec. HVI our experimen- 
tal results for the evolution of the I — V characteristics 
as a function of temperature and stack thickness. We 
compare results obtained from vertical transport through 
stacks of layers to in-plane measurements on monolayers 
and tetralayers and discuss these data within the context 
of recent models for charging-energy-mediated tunneling 
transport. Sec. IVl contains brief conclusions. 



II. CURRENT- VOLTAGE CHARACTERISTICS: 
BACKGROUND 

The purpose of this section is to discuss the general fea- 
tures of the current-voltage characteristics predicted for 
nanoparticle arrays in the semi-classical Coulomb block- 



ade regime and in the cotunneling regime. The semiclas- 
sical or " orthodox theory" of Coulomb blockade predicts 
exponential suppression of the conductance at low tem- 
peratures! 7 - In the zero-bias limit this suppression leads 
to go(T) = g exp[— Ec / ksT] for transport through a sin- 
gle nanoparticle connected by tunnel junctions to elec- 
trodes. Here E c is the charging energy, for a weakly cou- 
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pled particle approximated by Ec = a , where e is 
the dielectric constant, eo the permittivity of free space, 
and a is the particle radius. The bare conductance g 
corresponds to the temperature-independent Ohmic tun- 
nel conductance a junction would have in the absence of 
all charging effects. For nanoparticles, small a implies 
large Ec and thus a significant suppression of sequential 
tunneling at low temperatures. Conduction then occurs 
by one of two means: At sufficiently large applied bias 
voltage, the Coulomb blockade is overcome so that se- 
quential tunneling can commence. In arrays, this hap- 
pens for V > Vt, where Vt is a global Coulomb threshold 
voltage discussed below. Conversely, while a Coulomb 
blockade prevents direct sequential tunneling for V < Vt , 
cotunncling becomes a viable conduction channel. We 
note here that, for logic or switching devices based on 
single electron tunneling, cotunneling typically needs to 
be avoided. In principle, arrays are well-suited for this 
purpose since the probability for system-spanning cotun- 
nel events decreases dramatically with increasing system 
array size while Vt increases linearly . 17 ' 33 Nevertheless, as 
we show here, multiple shorter cotunnel events still can 
contribute significantly to the overall conduction. 



A. Semi-classical regime, V > Vt 

Applying the semiclassical picture to large arrays of 
nanoparticles with a random distribution of local chemi- 
cal potentials, Middleton and Wingreen showed that, at 
T = OK, conduction only occurs beyond a global thresh- 
old voltage, Vt = aNEc, where N is the number of 
nanoparticles spanning the gap between the electrodes, 
and a is a prefactor around 0.2-0.5 that depends on di- 
mensionality and array geometry^ 3 - At T — 0, in this 
model there is no conductance for V < Vt- Above Vt, the 
electrons percolate through the array via a multitude of 
branching paths that navigate local Coulomb blockade 
thresholds and optimize the total charging energy cost. 
Non-linear I — V characteristics with a power law depen- 
dence on V — Vt emerge as a consequence of the branch- 
ing^ 3 . 

The semiclassical T = scenario for transport above a 
global Coulomb threshold can be extended to finite tem- 
peratures^ As long as the charge disorder is sufficiently 
strong (an assumption easily satisfied in nanoparticle ar- 
rays because Ec is large) , the main effect of temperature 
is to wash out differences between chemical potentials 
of neighboring particles and thus " erase" local Coulomb 
blockade thresholds. As a consequence, the global thresh- 
old Vt decreases linearly with temperature, while the 
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power law exponent in the I — V characteristics remains 
unaffected. 

This picture is supported by experimental results on a 
variety of metal nanoparticle systems.— i 11 ' 12 Specifically 
for two-dimensional (2D) close-packed Au particle arrays 
it was found that V t (T) = V t (l - T/T*) where V t = 
aNEc with a = 0.23 and a temperature-independent 
exponent close to 2.2<&il The weak, linear dependence 
of the global threshold on temperature is very robust 
and persists in the presence of structural disorder in the 
array. On the other hand, such disorder leads to J — 
V characteristics that deviate from a simple power law 
above V t (T)&±2. 

By design, the Middleton and Wingren model and its 
extensions only apply to the large bias limit and do not 
consider thermally activated charge excitations. The low- 
bias conductances across "non-erased" local thresholds 
are assumed to be exponentially suppressed and, as a 
result, the global low-bias conductance is effectively con- 
trolled by the largest threshold along the optimal path. 
Because this path can avoid high charging energy costs to 
some extent by branching, for 2D close-packed arrays the 
relevant activation energy is not the full charging energy, 
but instead U « 0.2Ec and which leads to a zero-bias 
conductance go(T) = ,gexp[— U/kBT]A^ 7 - This simplified 
picture neglects processes that involve charge transfer 
over distances larger than a single particle (such as the 
cotunneling processes discussed next), but it is applica- 
ble at sufficiently high temperatures where the zero-bias 
conductance crosses over to simple activated, Arrhenius 
behavior characteristic of nearest neighbor hopping. 



B. Cotunneling regime, V < Vt 

Cotunneling, first discussed by Averin and Nazarov, 
is a tunneling process whereby electronic charge is trans- 
ferred through several neighbor particles cooperatively 2 ^.. 
For sequential tunneling each electron hop is a separate 
quantum event. By contrast, in cotunneling all charge 
transfers from the initial particle to the n th final particle 
occur via virtual states and thus count as a single quan- 
tum eventi. In principle, there are two different types 
of cotunneling processes: elastic and inelastic. In elastic 
cotunneling, an electron of the same eigenstate tunnels 
from the initial to the final particle, while inelastic cotun- 
neling involves the cooperative motion of multiple elec- 
trons, tunneling from one particle to the nearest neighbor 
in concert. Because elastic cotunneling becomes signifi- 
cant only at very low temperatures (below w IK in our 
system) we will focus in the following on inelastic pro- 
cesses. 

For the low-bias regime, the key quantity to con- 
sider for calculating the conductance is the probabil- 
ity to create an electron-hole excitation from an ini- 
tially neutral particle. This probability is proportional 
to exp[— E/ksT], where E is the electro-static energy as- 
sociated with the electron-hole pair. If charge transport 



occurs by sequential tunneling, the electron-hole pair is 
separated by one particle after the event and E is simply 
Ec- For cotunneling the final separation can be larger. 
If, at the end of an inelastic cotunneling process along 
a chain of neighboring particles, j electrons have each 
moved one particle down the current path, the electron- 
hole pair is separated by j + 1 particles, and E — As 
a result, for ksT « Ec cotunneling will always give rise 
to a larger zero-bias conductance than sequential tunnel- 
ing (even though cotunneling is a higher order process 
whose amplitude scales as (g/g q y , see below )2i2I. 

For the limit of zero applied bias, it was shown 2 -^ 
that cotunneling events span a typical distance, r* . In 
terms of the number of tunnel junctions involved, this 
distance is given by r* /d = (jr^ %) 1 ^ 2 , where £ is the 
localization length (£ < 2a for weakly coupled particles) 
and d the particle center-to-center spacing {d ks 8nm in 
our arrays). Equating this distance with j in the formula 
above for the activated creation of a dipole excitation 
leads to 



9o( T ) = 5exp(- v /T /T), 



(la) 



for the zero-bias conductance, with a characteristic tem- 
perature Tq given by 
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where C ~ 2. 8. 35 This result has the same functional 
form for the temperature dependence as Efros-Shklovskii 
variable range hopping (VRH) for doped semiconduc- 
tors.— It similarly describes the effect of balancing net 
tunnel distance and energy cost. However, it arises here 
not from direct, single charge tunnel events between 
sparsely distributed atomic defect sites, but from coop- 
erative multi-electron processes in dense arrays of metal 
particles. 

At finite but small voltage, for inelastic cotunneling 
processes involving a distance of j junctions, the current 
is given by 
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Here Vj c t is the bias voltage drop across each of the j 
junctions, g and g q are the junction and quantum con- 
ductances respectively, and E — the energy cost for 
the dipole creation. For Vj Ct << k^T the conduction 
is Ohmic and optimization with respect to j returns the 
result for go(T) in Eq. (fTa)) . 

If the current path involves sufficiently many particles 
in series, we expect that multiple cotunnel events will oc- 
cur, each over a typical distance j = r* jd. This distance 
will increase according to T -1 / 2 as the temperature is 
reduced, leading to a corresponding increase in the cur- 
vature of the low-bias I — V characteristics. On the other 
hand, if the length of the current path is fixed at a dis- 
tance shorter than r*/d, the temperature dependence of 
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the curvature will saturate. This paper explores these 
effects by performing transport measurements in metal 
nanoparticle arrays of different thicknesses in the direc- 
tion perpendicular to the array plane. We compare these 
results with those of monolayers and tetralayers where 
current-voltage measurements were taken in-plane. 



III. EXPERIMENTAL PROCEDURE 

Dodecanethiol-ligated gold nanoparticles were synthe- 
sized by the digestive ripening method described in 
Refi^S. This provided particle diameters around 5.5nm 
with tight size control and dispersion less than 5% in 
each batch. The nanoparticle arrays were deposited 
onto silicon substrates coated with lOOnm amorphous 
Si^N^. For sample characterization by transmission 
electron microscopy (TEM), these substrates contained 
70/im x70/im or 300/im x 300/im "window" areas under 
which the Si had been etched away to leave free-standing, 
TEM-transparent S13N4 membranes. Before particle de- 
position, we created by electron beam lithography and 
thermal evaporation 20nm thick chromium electrodes on 
top of the substrates, reaching into the window areas. For 
in-plane measurements, sets of two opposing electrodes 
were fabricated, 2/xm in width and with gap 500nm be- 
tween them. For transport measurements perpendicular 
to the array plane, the bottom electrode consisted of a 
Cr strip 10/im in width and 600/im in length. After par- 
ticle deposition, a top electrode was evaporated over the 
array, consisting of a Ni strip 30nm in thickness, 300/im 
in length and 10/im in width. The top electrodes were 
shadow evaporated through masks made from the same 
Si%N4 window substrates used for the samples, with CF4 
etched rectangular holes of dimension 300/im x 10/im. 
By orienting the top electrodes orthogonal to the bottom 
electrodes in a cross-like geometry, overlapping regions of 
area 10/im x 10/im were created. Cr and Ni were chosen 
as the electrode material because of their excellent ad- 
hesion to the SizN± as well as to the gold particles, and 
because uniform, continuous electrodes can be formed 
reliably at relatively low thickness. The use of Ni kept 
the heat load on the array during evaporation lower than 
would have been the case with Cr. However, the upper- 
most dodecancthiol ligands were likely burnt away during 
this process, bringing the top electrode in direct contact 
with top particle layer. This is substantiated by the ver- 
tical transport measurements discussed below. On the 
other hand, we were able to ascertain by TEM that the 
top electrode evaporation did not lead to any noticeable 
sintering of the particles. Figure [1] displays an optical im- 
age of a bilayer sample with top and bottom electrodes. 

Gold nanoparticle monolayers were created at the air 
interface of a water droplet, similar to a Langmuir- 
Schaefer technique 3 *^. First, a ~ 5mm diameter water 
droplet was created on a glass slide. This was followed 
by deposition of 20-30/iL of gold nanoparticle solution 
(concentration 10 13 particles/mL) on top of the droplet. 




FIG. 1: Optical image of a bilayer nanoparticle array on a 
S13N4 window substrate. The top and bottom electrodes 
(light rectangular strips) are oriented perpendicular to one 
another. The large dark rectangle is the window area. 



Using this technique, the particles quickly spread to form 
a beautiful, compact monolayer across the water surface, 
presumably due to surface tension. By dipping an in- 
verted substrate with prefabricated electrodes onto the 
top of the droplet, the monolayer was then transfer- 
printed onto the substrate. This preparation method 
gave interparticle distances of 1.7nm to 1.8nm between 
different batches of nanoparticle solutions. Multilayer ar- 
rays were created via layer- by- layer deposition, repeating 
the process described above. There is no registry between 
subsequent layers using this technique, but the compact, 
highly ordered particle arrangement is maintained. Fig- 
ure [2^ shows a monolayer spanning the region between 
two planar electrodes. TEM images of multilayer arrays 
for vertical transport measurements in Figs. [2)>d were 
taken right next to the electrodes because the top elec- 
trodes in these samples were not TEM-transparent. The 
data presented in this paper was taken on mono-, bi-, tri, 
and tetralayers. 

Current- voltage (/ — V) measurements were performed 
with a Keithley 6430 source meter. Bias voltages of V < 
\4V\ were applied perpendicular to the sample plane at 
rates of 0.05-1.5mV/sec and at temperatures 5K-100K. 
For in-plane I—V measurements, bias voltages V < \20V\ 
were applied at rates of 5-25mV/sec over the range from 
10K - 160K. The cryostat was wired with coax cables 
for low-leakage, low-noise measurements. Typical current 
noise was below 40fA. TEM images of all samples were 
taken posterior to I — V measurements to preclude any 
irradiation damage or contamination of the samples, even 
if undetectable in the image, from affecting the transport 
measurements . 

The zero-bias conductance was obtained by perform- 
ing I — V measurements, at fixed temperature, covering 
a narrow voltage range across the origin and extracting 
go(T) from the slope. Care was taken to restrict this volt- 
age range to the linear response. Because of the current 
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FIG. 2: Transmission electron microscopy (TEM) images of 
nanoparticle arrays used in the experiments: (a) monolayer, 
(b) bilayer, (c) trilayer, and (d) tetralayer. The main pan- 
els show close-up views of the highly ordered particle ar- 
rangement. The bottom insets in a-d display diffraction pat- 
terns computed by fast Fourier transform. The top insets are 
zoomed-out images of the area covering the region between 
the in-plane electrodes of the monolayer (a), and the region 
near the bottom electrodes of bi-, tri-, and tetra-layers mea- 
sured in perpendicular direction (b-d). The dark regions of 
in these insets are the electrodes. 

noise floor, this meant that go(T) could be tracked reli- 
ably down to temperatures around 20K, but not always 
below that value. This same issue also restricted the ex- 
tent over which we could track the nonlinear behavior at 
the lowest temperatures before becoming limited by the 
noise floor. 

Transport measurements in the presence of the mag- 
netic field were performed in a 14/16T Oxford He-3 cryo- 
stat. Magnetic fields up to 10T were applied both paral- 
lel and transverse to the sample plane at temperatures of 
5K, IK and 400mK. The voltage sweep rate was 0.5 — 1.2 
mV/s. 



IV. RESULTS AND DISCUSSION 

A. Non-Ohmic Regime 

We first focus on the regime of finite applied bias. Here 
the transport properties are highly non-Ohmic. Figs. [3] 
and[l]show our I — V data for vertical and in-plane trans- 
port, respectively. In both figures, the left and right 
columns show the same data, but on linear and on double 
logarithmic scales. By appropriately normalizing the cur- 
rent so that I — V traces taken at different temperatures 
can be visualized on the same linear plot, the increase 
in nonlinearity with decreasing temperature is immedi- 



ately apparent. This is the purpose of the left columns of 
panels. The right columns show the raw data and indi- 
cate how the highly nonlinear and strongly temperature- 
dependent behavior at finite, but small, bias merges into 
a much less temperature-dependent high bias regime. 
For in-plane transport across large arrays this high bias 
regime corresponds to the semiclassical V > Vt regime. 
Here, all / — V traces are predicted to approach the same 
power law behavior with exponent close to 2, merging 
onto this asymptote at different temperatures given by 
V t (T). This is borne out by the data in Fig. gj;,d. For 
vertical transport, the distance traversed at high bias is 
the thickness of the stack of layers (a few particle dis- 
tances) , and this is too short for a path branching pattern 
to emerge. Instead, the behavior seen in Fig. [3] resembles 
that of short parallel paths, each retaining the strongly 
temperature-dependent prefactor in Eqj2] (as opposed to 
acquiring the weak, linear T-dependence of Vt(T)). 



1. Vertical Transport 

As can be seen from data on bi-, tri- and tetra- 
layers (Fig. [3^i-c), the curvature of the small bias, low- 
temperature I — V characteristics increases with increas- 
ing thickness of the stack. The data are fit well by power 
laws of the form / ~ V 13 , with exponents f3 as indicated 
in the figure. At 5K, the exponents are approximately 3, 
5, and 7 for a bi-, tri- and tetralayer, respectively. While 
(3 for the bilayer remains the same as the temperature 
increases to 10K, /3 for tri- and tetralayers decreases to 4 
and 5. 

In order to compare these data with theoretical pre- 
dictions, we first estimate the charging energy. Using 
Ec = 47ree - , particle radius a ss 2.5nm as measured by 
TEM and e » 4 we find E c « 125meV or 1600K. Here 
the value of the dielectric constant is based on our pre- 
vious results on multilayers.— 

For vertical transport across a short stack of n parti- 
cles, the applied bias voltage and the bias per junction are 
related by Vj c t = — • Based on Eq. @ we therefore expect 
that the nonlinear cotunneling regime ksT < eVj c t < Eq 
should persist up to a good fraction of a Volt applied 
bias at low temperatures. This is indeed what we find in 
Fig. [3J Equating the exponents (3 found from the power 
law fits to the exponent 2j — 1 in Eq. © we find the 
typical number of junctions participating in the cotunnel 
events. At the lowest temperatures, this leads to 2, 3 and 
4 junctions for bi-, tri- and tetralayers, respectively, con- 
sistent with the notion that the top electrode has fused 
directly to the uppermost particle layer. The fact that (3 
for the bilayer saturates at a value around 3 below 10K is 
an indication that the stack thickness provides a cut-off: 
dipole separations larger than 2 junctions are not possi- 
ble because they would exceed the distance between the 
top and bottom electrodes (however, see the caveat for 
the zero-bias limit, below). Consequently, at these low 
temperatures all cotunneling events span the electrode 
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FIG. 3: Current-voltage (I — V) characteristics as a function 
of temperature for bilayer (a), trilayer (b), and tetralayer (c) 
arrays measured in the direction perpendicular to the array 
plane. Panels d, e and f are log-log plots of the data in pan- 
els a, b and c, respectively; the solid lines are guides to the 
eye corresponding to power laws with exponents as indicated. 
The current scales in panels a-c are normalized with respect 
to the current at fixed bias voltage (0.4V in (a), 0.6V in (b), 
0.9V in (c)) to bring out the increasing nonlinearity as tem- 
perature is reduced. The solid lines in panels a-c correspond 
to power law fitting curves of the form I ~ V B . The exponents 
[3 obtained from these fits have uncertainties of ±0.1to ± 0.2 
in panels a and b, and ±0.3 in panel c. 



distance, the exponential term accounting for the ther- 
mal creation of the electron-hole pair inside the array 
no longer remains, and Eq. ([2|) reduces to the the form 
known from short linear arrays of mesoscopic isfand o 28 ! 32 
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where the term in the last square brackets was 
taincd by expanding [1 ± (^4^-) 2 p _1 . As before, 
power-law exponent equals 2j — 1 in the voltage range 
kgT < £|- < Ec and reduces to 1 at high temperatures. 

Using Eq. [3J appropriate for small , we find that the 
current between 4.2K and 10K is predicted to increase by 
a factor of 6. This is in reasonable agreement with the 
data which shows a factor of 8-10. For the tri- and the 
tctralayers, the cotunneling length is on the order of the 
sample length only at 5K. At higher temperatures, this 
length becomes less than the distance between the elec- 
trodes, and thus multiple cotunneling will be involved. 




2.0 r 

1.5- 
1.0 

0.5- 
0.0- 





/ b 


I -V s " 




" |-v" 


I 


- I~V 33 






J° 10K 


l-V" 


Jl ° 30K 


- I-V" 


» 55K 




o 85K 




* 120K 


1 2 


3 4 5 6 7 











2 




.... ss5p 




1 0K 




o 


30K " 






50K 




O 


80K ' 




. . . ° 


120K 


10° 


10 1 




V(V) 










— j- 




2 






□ 


10K 




o 


30K 




A 


55K ; 




o 


85K 




-* 


120K . 



10° 
10" 

3 

10 2 r 

10° 
10' a 

10 6 
10 4 

10' I 

10° 

10 2 



V(V) 



V(V) 



FIG. 4: Current-voltage (I — V) characteristics as a function 
of temperature for monolayer (a) and tetralayer (b) arrays 
measured along the array plane. Panels c, d are log-log plots 
of the data in panels a, b, respectively. The current scales in 
panels a and b are normalized with respect to the current at 
fixed bias voltage of 5V to bring out the increasing nonlinear- 
ity as temperature is reduced. The solid lines in panels a, b 
correspond to power law fitting curves of the form I ~ V 13 . 
The exponents /3 obtained from these fits have uncertainties 
of ± 0.2 to ±0.3. 



2. In-plane Transport 

I — V measurements taken in-plane for comparison ex- 
hibit behavior similar to the perpendicular measurements 
just discussed. This is shown in Figs. [5] for a mono- and a 
tetralayer (for additional data see Refj22). Because these 
are large arrays where the number of junctions across 
the gap is ~ 70, the global Coulomb blockade threshold 
Vt ~ NEq easily reaches several Volts. As for vertical 
transport, the data in the small-bias regime below Vt are 
fit well by power laws of the form / ~ V@, with exponents 
(3 as indicated in the figure. Equating (3 with the expo- 
nent 2j — 1 in Eq. ^ similarly gives the typical cotun- 
nel distance j as a function of temperature. The power- 
law exponents for both the monolayers and tetralayers 
measured in-plane imply that the maximum number of 
particles participating in one cotunneling event is much 
smaller than the number of particles spanning the gap 
between the electrodes, thus necessitating multiple co- 
tunnel events to get across. The fact that the tetralayer's 
power-law exponent (3 ~ 5 at 10-ftT is similar to the value 
for the same exponent measured at 10K perpendicular 
to the array plane (~ 4.2), supports the notion that the 
same cotunnel mechanism is at work leading to a similar 
typical cotunnel distance in both cases. For V > Vt, all 
I — V curves turn over to approach a power-law exponent 
close to 2 as predicted by the semiclassical picture, and 
for V approaching the zero-bias limit, the I — V traces 
become linear, in line with Eq. 
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B. Ohmic Regime 

Similar to our observations on samples measured in- 
plane^, go(T) in the perpendicular direction also fol- 
lows the stretched exponential form given by Eqs. (fla|) 
and (|lb[) . As mentioned earlier, this behavior echoes 
VRH, but here the same functional form with the char- 
acteristic inverse square root temperature dependence in 
the exponential arises from multiple cotunneling. Inter- 
estingly, we find this functional form for all stack thick- 
nesses. In interpreting this result, one potential problem 
is that, for the vertical transport case discussed here, a 
straight path form the bottom to the top electrode in- 
volves only a few junctions (no more than perhaps 4 for 
the tetralayer). Therefore, at first glance multiple co- 
tunnel events and the optimization of cotunnel distance 
versus energy cost, both necessary to produce the inverse 
square root temperature dependence of go (T) , appear un- 
likely. However, it is reasonable that the paths taken by 
the electrons depend on the presence of the applied field. 
As V — ► 0, the electrons might be able to meander into 
the lateral, in-plane direction to find an optimum, ener- 
getically cost-effective path to tunnel from the bottom to 
the top electrode. At higher applied bias, the electric field 
channels the electrons along a more direct path, and thus 
the cotunnel distance cannot exceed the number of tunnel 
junctions across the stack thickness. We speculate that 
this scenario is the reason why the low-temperature I — V 
characteristics, measured at finite bias, indicate a rela- 
tively small number of junctions commensurate with the 
stack thickness, while the zero-bias conductance, mea- 
sured at a significantly smaller bias (feasible at higher 
temperatures), hints at much longer paths. 

From the slopes of the semi-log plots in Fig. [5l we find 
the characteristic temperatures, To and, using Eq. Ilbl 
the localization lengths, £ for bi-, tri- and tetra-layers. 
Table U compares To values from both in-plane and per- 
pendicular measurements (including data from Ref. 1321). 
All values are between 72-83K, except for vertical trans- 
port through a bilayer. The T in this case is significantly 
smaller, implying correspondingly larger values for ei- 
ther the localization length or the dielectric constant (see 
Eg. Ilb[) . Since the localization lengths obtained from To 
for all other systems are w 2nm, and thus less than a 
particle diameter as appropriate for cotunneling through 
weakly coupled particles, it seems reasonable that the re- 
duced T) is associated with an increase in the effective e. 
This is likely to be produced by the close proximity of the 
top electrode required for vertical transport. As soon as 
the electrodes are 1 to 2 layers away, the dielectric con- 
stant of the systems seems to be unaffected, i.e., the To 
of both trilayers and tetralayers are similar. The same 
argument might also explain why the in-plane monolayer 
exhibits a slightly larger To than most other systems: it 
experiences less screening from additional layers above it. 

Using r/d = (-^jt f )^ 2 together with Eq. (fTb|) we can 
estimate the typical cotunnel distances from the exper- 
imental values for Tq. In Fig. [S] we plot the result in 




FIG. 5: Temperature-dependent zero-bias conductance of a 
monolayer measured in-plane, and of a bilayer, trilayer and 
tetralayer where I — V measurements were orthogonal to the 
sample plane. The inset shows Arrhenius behavior in the 
monolayer with I — V measurements performed in the lat- 
eral direction, and in the trilayer and tetralayer with I — V 
measurements taken perpendicular to the sample plane. 



Sample 


T (K) (in- plane) 


To{K) (Perpendicular) 


Monolayer 


6.5 x 10 3 




Bilayer 


3.3 x 10 3 


2 x 10 3 


Trilayer 


5.1 x 10 3 


6.4 x 10 3 


Tetralayer 


5.2 x 10 3 


5.8 x 10 3 



TABLE I: Comparison of the characteristic temperatures To 
obtained from fits of the zero-bias conductance data to Eg. Ilbl 
To values have uncertainties of ~ 5%. 



terms of the number of junctions involved, j — r*/d — 
(ct0 1/2 §r (dotted lines), where C w 2.8. In the same 
figure, this inverse-square-root temperature dependence 
obtained from the zero-bias conductance is compared to 
the j-values extracted directly from the curvature of the 
I — V characteristics at finite bias. The two sets of data 
track each other well, differing by less than one junction 
over most of the measured temperature range. For exam- 
ple, at 10K, the distances determined based on To are 4 
junctions for both tri- and tetra-layers, while the power- 
law exponents give 2.5 and 3 junctions, respectively. 

At sufficiently high temperatures, the cotunnel dis- 
tance approaches a single junction. At this point go(T) 
crosses over to simple activated, Arrhenius behavior, with 
5o (T) = gexp[— -^f] (dotted lines in Fig. [5] and plot in 
inset). For vertical transport in both tri- and tetra-layers 
this occurred above roughly 60K and led to U/ks ~ 
193K and 228K, but in the bilayer was not observed over 
the full measured range up to 125K. For comparison, in- 
plane transport data for the monolayer shown in Fig. [5] 
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FIG. 6: Typical cotunnel distance as a function of temper- 
ature for trilayer (a) and tetralayer (b) arrays measured in 
perpendicular direction, and for a monolayer (c) measured 
in-plane. The distance is given as the number j of junctions 
along a chain of particles involved in the cooperative cotunnel 
process. The plots compare the dependence extracted from 
the zero-bias conductance (dotted lines) with the values for j 
obtained directly from the I — V power-law exponents in the 
small bias regime. 



gives a cross-over temperature ~ 70K and ll/kg ~ 300K. 
Similar cross-overs at temperatures around 70-100K were 
observed in our previous previous data for the in-plane 
zero-bias conductance of bi-, tri, and tetralayers^ From 
the expressior&iL 3 - 2 U ~ 0.2NEc that relates the mea- 
sured high-temperature, Arrhenius-type activation en- 
ergy U to the charging energy we find Ec/ks values of 
965K, 1140K and 1500K for the tri-, tetra-, and mono- 
layer in Fig. [5l respectively. This is in reasonable agree- 
ment with the value of 1600K we estimated above using 

2 

Ec = a ' es P ecian y given the experimental uncer- 
tainty in determining U over the available limited tem- 
perature range. 



C. Magneto-transport 

In doped semiconductors, variable-range hopping at 
low bias voltage results in a change of conductance in 
the presence of a magnetic field. Large magnetic fields 
squeeze the exponentially-decaying wavefunctions in the 
transverse direction, i.e. if the applied magnetic field 
is along the out-of-plane or z-direction, the wavefunc- 
tions are squeezed in the x- and y-directions. The de- 
crease in the wavefunction overlap gives rise to a de- 
creased tunneling amplitude, and this in turn produces 
a decrease in the variable-range- hopping conductance 35 . 



FIG. 7: Low-bias magnetotransport data for a monolayer 
measured in-plane at 400mK (a) and for a bilayer measured in 
perpendicular direction at 5K (b). In both cases the field was 
applied perpendicular to the sample plane, i.e., perpendicular 
to the current flow in (a) and parallel to it in (b). 



This phenomenon has been observed in both n-type 36 
and p-type semiconductors^. Recent magneto-transport 
data on n-type CdSe semiconductor nanoparticle films 
has also been interpreted within the VRH picture^!. By 
investigating the magnetic field dependence of our metal 
nanoparticle arrays, we therefore can provide a further 
test for the transport mechanism: VRH is expected to 
produce a positive magnetoresistance, while inelastic co- 
tunneling should not be affected even by large magnetic 
fields since the hops of all participating electrons involve 
only nearest neighbors. 

Figure [7] shows that the presence of a magnetic field 
perpendicular to the sample plane has no clear measur- 
able effect within our experimental resolution. These 
data were taken for in-plane transport through a mono- 
layer and for vertical transport through a bilayer. The 
same behavior was observed when magnetic fields were 
applied along the direction of the sample plane. This 
magnetic-field-independent of the I — V characteristics 
is to be compared to the large positive magnetoresis- 
tance (~ 150% at I0T and T = 0.3K) reported for 
CdSe films. 3 — The distance over which an applied mag- 
netic field H affects charge motion is set by the mag- 
netic length L = ($ /-£0 5 , where $ = 55 is the flux 
quantum. Since L > 14nm at fields H < I0T the ab- 
sence of magnetoresistance in the weakly coupled metal 
nanoparticle arrays discussed here implies that the mo- 
tion of individual charges occurs over distances less than 
L, in line with inelastic cotunneling. At sufficiently low 
temperatures where a cross-over to elastic cotunneling 
takes place, magnetoresistance should reappear. In our 
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samples, this cross-over is expected below IK, and it 
may be possible that this produced the slight differences 
between the traces with and without magnetic field in 
Fig. [T^i. However, we believe that the data are too close 
to the experimental noise floor. 

V. CONCLUSION 

The results presented here demonstrate that the elec- 
tronic transport of weakly-coupled metal nanoparticle ar- 
rays is well described by a picture based on sequential 
tunneling in the high-bias regime above a global Coulomb 
blockade threshold, and by inelastic cotunncling in the 
low-bias regime below that threshold. In particular, mul- 
tiple inelastic cotunneling provides an explanation for 
both the exponential, inverse-square-root temperature 
dependence of the conductance, go(T), in the limit of 
zero bias and the non-linear, power law I — V charac- 
teristics observed inside the Coulomb-blockade regime 
at finite bias. While the zero-bias temperature depen- 
dence of go(T) mimics that known from single-charge, 
variable-range hopping transport in doped semiconduc- 
tors, cotunneling events are cooperative, multi-electron 
processes. In close-packed metal particle arrays individ- 
ual electrons are unlikely to move further than to their 



nearest neighbor, at least at temperatures above IK, in 
contrast to hopping between sparse doping sites. Com- 
parison of transport measurements taken in-plane and 
perpendicular to the arrays (Figs. [3] and Ql indicate that 
typical cotunnel events in Au nanoparticle arrays involve 
up to 4 electrons (and thus reach net distances corre- 
sponding to 4 tunnel junctions in a row) at 10K (Fig. [5]). 
The absence of clean evidence for magnetoresistive effects 
up to the highest applied magnetic fields (10T) further 
supports the picture of inelastic cotunneling as the mech- 
anism for charge transport inside the Coulomb blockade 
regime. 
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